Dissection of TAF1 neuronal splicing and implications for neurodegeneration in X-linked dystonia-parkinsonism.

X-linked dystonia-parkinsonism (XDP) is a monogenic neurodegenerative disorder of the basal ganglia, which presents as a combination of hyperkinetic movements and parkinsonian features. The underlying genetic mechanism involves the insertion of a SINE-VNTR-Alu retrotransposon within the TAF1 gene. Interestingly, alterations of TAF1 have been involved in multiple neurological diseases. In XDP, the SINE-VNTR-Alu insertion in TAF1 has been proposed to result in alternative splicing defects, including the decreased incorporation of a neuron-specific microexon annotated as 34'. This mechanism has become controversial as recent studies failed to provide support. In order to resolve this conundrum, we examined the alternative splicing patterns of TAF1 mRNAs in XDP and control brains. The impact of the disease-associated SINE-VNTR-Alu on alternative splicing of microexon 34' was further investigated in cellular assays. Subsequently, microexon 34' incorporation was explored by RT-PCR and Nanopore long-read sequencing of TAF1 mRNAs from XDP and control brains tissues. Using cell-based splicing assays, we demonstrate that presence of the disease-associated SINE-VNTR-Alu does not affect the inclusion of microexon 34'. In addition, we show that (1) microexon 34'-containing TAF1 mRNAs are detected at similar levels in XDP as in controls and that (2) the architecture of TAF1 transcripts is remarkably similar between XDP and controls brains. These results indicate that microexon 34' incorporation into TAF1 mRNA is not affected in XDP brains. Our findings shift the current paradigm of XDP by discounting alternative splicing of TAF1 microexon 34' as the molecular basis for this disease.

Evolution of regulatory signatures in primate cortical neurons at cell-type resolution.

The human cerebral cortex contains many cell types that likely underwent independent functional changes during evolution. However, cell-type-specific regulatory landscapes in the cortex remain largely unexplored. Here we report epigenomic and transcriptomic analyses of the two main cortical neuronal subtypes, glutamatergic projection neurons and GABAergic interneurons, in human, chimpanzee, and rhesus macaque. Using genome-wide profiling of the H3K27ac histone modification, we identify neuron-subtype-specific regulatory elements that previously went undetected in bulk brain tissue samples. Human-specific regulatory changes are uncovered in multiple genes, including those associated with language, autism spectrum disorder, and drug addiction. We observe preferential evolutionary divergence in neuron subtype-specific regulatory elements and show that a substantial fraction of pan-neuronal regulatory elements undergoes subtype-specific evolutionary changes. This study sheds light on the interplay between regulatory evolution and cell-type-dependent gene-expression programs, and provides a resource for further exploration of human brain evolution and function.

Recently Evolved Enhancers Emerge with High Interindividual Variability and Less Frequently Associate with Disease.

Mutations in non-coding regulatory DNA such as enhancers underlie a wide variety of diseases including developmental disorders and cancer. As enhancers rapidly evolve, understanding their function and configuration in non-human disease models can have important clinical applications. Here, we analyze enhancer configurations in tissues isolated from the common marmoset, a widely used primate model for human disease. Integrating these data with human and mouse data, we find that enhancers containing trait-associated variants are preferentially conserved. In contrast, most human-specific enhancers are highly variable between individuals, with a subset failing to contact promoters. These are located further away from genes and more often reside in inactive B-compartments. Our data show that enhancers typically emerge as instable elements with minimal biological impact prior to their integration in a transcriptional program. Furthermore, our data provide insight into which trait variations in enhancers can be faithfully modeled using the common marmoset.

Hominin-specific regulatory elements selectively emerged in oligodendrocytes and are disrupted in autism patients.

Speciation is associated with substantial rewiring of the regulatory circuitry underlying the expression of genes. Determining which changes are relevant and underlie the emergence of the human brain or its unique susceptibility to neural disease has been challenging. Here we annotate changes to gene regulatory elements (GREs) at cell type resolution in the brains of multiple primate species spanning most of primate evolution. We identify a unique set of regulatory elements that emerged in hominins prior to the separation of humans and chimpanzees. We demonstrate that these hominin gains perferentially affect oligodendrocyte function postnatally and are preferentially affected in the brains of autism patients. This preference is also observed for human-specific GREs suggesting this system is under continued selective pressure. Our data provide a roadmap of regulatory rewiring across primate evolution providing insight into the genomic changes that underlie the emergence of the brain and its susceptibility to neural disease.

Neuronal-specific microexon splicing of TAF1 mRNA is directly regulated by SRRM4/nSR100.

Neuronal microexons represent the most highly conserved class of alternative splicing events and their timed expression shapes neuronal biology, including neuronal commitment and differentiation. The six-nt microexon 34' is included in the neuronal form of TAF1 mRNA, which encodes the largest subunit of the basal transcription factor TFIID. In this study, we investigate the tissue distribution of TAF1-34' mRNA and protein and the mechanism responsible for its neuronal-specific splicing. Using isoform-specific RNA probes and antibodies, we observe that canonical TAF1 and TAF1-34' have different distributions in the brain, which distinguish proliferating from post-mitotic neurons. Knockdown and ectopic expression experiments demonstrate that the neuronal-specific splicing factor SRRM4/nSR100 promotes the inclusion of microexon 34' into TAF1 mRNA, through the recognition of UGC sequences in the poly-pyrimidine tract upstream of the regulated microexon. These results show that SRRM4 regulates temporal and spatial expression of alternative TAF1 mRNAs to generate a neuronal-specific TFIID complex.

Transcriptomic and Epigenomic Profiling of Histone Deacetylase Inhibitor Treatment Reveals Distinct Gene Regulation Profiles Leading to Impaired Neutrophil Development.

The clinical use of histone deacetylase inhibitors (HDACi) for the treatment of bone marrow failure and hematopoietic malignancies has increased dramatically over the last decades. Nonetheless, their effects on normal myelopoiesis remain poorly evaluated. Here, we treated cord blood derived CD34+ progenitor cells with two chemically distinct HDACi inhibitors MS-275 or SAHA and analyzed their effects on the transcriptome (RNA-seq), epigenome (H3K27ac ChIP-seq) and functional and morphological characteristics during neutrophil development. MS-275 (entinostat) selectively inhibits class I HDACs, with a preference for HDAC1, while SAHA (vorinostat) is a non-selective class I/II HDACi. Treatment with individual HDACi resulted in both overlapping and distinct effects on both transcriptome and epigenome, whereas functional effects were relatively similar. Both HDACi resulted in reduced expansion and increased apoptosis in neutrophil progenitor cells. Morphologically, HDACi disrupted normal neutrophil differentiation what was illustrated by decreased percentages of mature neutrophils. In addition, while SAHA treatment clearly showed a block at the promyelocytic stage, MS-275 treatment was characterized by dysplastic features and skewing towards the monocytic lineage. These effects could be mimicked using shRNA-mediated knockdown of HDAC1. Taken together, our data provide novel insights into the effects of HDAC inhibition on normal hematopoietic cells during neutrophil differentiation. These findings should be taken into account when considering the clinical use of MS-275 and SAHA, and can be potentially utilized to tailor more specific, hematopoietic-directed HDACi in the future.

AML Subtype Is a Major Determinant of the Association between Prognostic Gene Expression Signatures and Their Clinical Significance.

Relapse in acute myeloid leukemia (AML) may result from variable genetic origins or convergence on common biological processes. Exploiting the specificity and sensitivity of regulatory DNA, we analyze patient samples of multiple clinical outcomes covering various AML molecular subtypes. We uncover regulatory variation among patients translating into a transcriptional signature that predicts relapse risk. In addition, we find clusters of coexpressed genes within this signature selectively link to relapse risk in distinct patient subgroups defined by molecular subtype or AML maturation. Analyzing these gene clusters and the AML subtypes separately enhances their prognostic value substantially and provides insight in the mechanisms underlying relapse risk across the distinct patient subgroups. We propose that prognostic gene expression signatures in AML are valid only within patient subgroups and do not transcend these subgroups.

Mapping the Global Chromatin Connectivity Network for Sox2 Function in Neural Stem Cell Maintenance.

The SOX2 transcription factor is critical for neural stem cell (NSC) maintenance and brain development. Through chromatin immunoprecipitation (ChIP) and chromatin interaction analysis (ChIA-PET), we determined genome-wide SOX2-bound regions and Pol II-mediated long-range chromatin interactions in brain-derived NSCs. SOX2-bound DNA was highly enriched in distal chromatin regions interacting with promoters and carrying epigenetic enhancer marks. Sox2 deletion caused widespread reduction of Pol II-mediated long-range interactions and decreased gene expression. Genes showing reduced expression in Sox2-deleted cells were significantly enriched in interactions between promoters and SOX2-bound distal enhancers. Expression of one such gene, Suppressor of Cytokine Signaling 3 (Socs3), rescued the self-renewal defect of Sox2-ablated NSCs. Our work identifies SOX2 as a major regulator of gene expression through connections to the enhancer network in NSCs. Through the definition of such a connectivity network, our study shows the way to the identification of genes and enhancers involved in NSC maintenance and neurodevelopmental disorders.

Intracranial Aneurysm-Associated Single-Nucleotide Polymorphisms Alter Regulatory DNA in the Human Circle of Willis.

BACKGROUND AND PURPOSE: Genome-wide association studies significantly link intracranial aneurysm (IA) to single-nucleotide polymorphisms (SNPs) in 6 genomic loci. To gain insight into the relevance of these IA-associated SNPs, we aimed to identify regulatory regions and analyze overall gene expression in the human circle of Willis (CoW), on which these aneurysms develop. METHODS: We performed chromatin immunoprecipitation and sequencing for histone modifications H3K4me1 and H3K27ac to identify regulatory regions, including distal enhancers and active promoters, in postmortem specimens of the human CoW. These experiments were complemented with RNA sequencing on the same specimens. We determined whether these regulatory regions overlap with IA-associated SNPs, using computational methods. By combining our results with publicly available data, we investigated the effect of IA-associated SNPs on the newly identified regulatory regions and linked them to potential target genes. RESULTS: We find that IA-associated SNPs are significantly enriched in CoW regulatory regions. Some of the IA-associated SNPs that overlap with a regulatory region are likely to alter transcription factor binding, and in proximity to these regulatory regions are 102 genes that are expressed in the CoW. In addition, gene expression in the CoW is enriched for genes related to cell adhesion and the extracellular matrix. CONCLUSIONS: CoW regulatory regions link IA-associated SNPs to genes with a potential role in the development of IAs. Our data refine previous predictions on SNPs associated with IA and provide a substantial resource from which candidates for follow-up studies can be prioritized.

Tomo-Seq Identifies SOX9 as a Key Regulator of Cardiac Fibrosis During Ischemic Injury.

BACKGROUND: Cardiac ischemic injury induces a pathological remodeling response, which can ultimately lead to heart failure. Detailed mechanistic insights into molecular signaling pathways relevant for different aspects of cardiac remodeling will support the identification of novel therapeutic targets. METHODS: Although genome-wide transcriptome analysis on diseased tissues has greatly advanced our understanding of the regulatory networks that drive pathological changes in the heart, this approach has been disadvantaged by the fact that the signals are derived from tissue homogenates. Here we used tomo-seq to obtain a genome-wide gene expression signature with high spatial resolution spanning from the infarcted area to the remote to identify new regulators of cardiac remodeling. Cardiac tissue samples from patients suffering from ischemic heart disease were used to validate our findings. RESULTS: Tracing transcriptional differences with a high spatial resolution across the infarcted heart enabled us to identify gene clusters that share a comparable expression profile. The spatial distribution patterns indicated a separation of expressional changes for genes involved in specific aspects of cardiac remodeling, such as fibrosis, cardiomyocyte hypertrophy, and calcium handling (Col1a2, Nppa, and Serca2). Subsequent correlation analysis allowed for the identification of novel factors that share a comparable transcriptional regulation pattern across the infarcted tissue. The strong correlation between the expression levels of these known marker genes and the expression of the coregulated genes could be confirmed in human ischemic cardiac tissue samples. Follow-up analysis identified SOX9 as common transcriptional regulator of a large portion of the fibrosis-related genes that become activated under conditions of ischemic injury. Lineage-tracing experiments indicated that the majority of COL1-positive fibroblasts stem from a pool of SOX9-expressing cells, and in vivo loss of Sox9 blunted the cardiac fibrotic response on ischemic injury. The colocalization between SOX9 and COL1 could also be confirmed in patients suffering from ischemic heart disease. CONCLUSIONS: Based on the exact local expression cues, tomo-seq can serve to reveal novel genes and key transcription factors involved in specific aspects of cardiac remodeling. Using tomo-seq, we were able to unveil the unknown relevance of SOX9 as a key regulator of cardiac fibrosis, pointing to SOX9 as a potential therapeutic target for cardiac fibrosis.

Circadian networks in human embryonic stem cell-derived cardiomyocytes.

Cell-autonomous circadian oscillations strongly influence tissue physiology and pathophysiology of peripheral organs including the heart, in which the circadian clock is known to determine cardiac metabolism and the outcome of for instance ischemic stress. Human pluripotent stem cells represent a powerful tool to study developmental processes in vitro, but the extent to which human embryonic stem (ES) cell-derived cardiomyocytes establish circadian rhythmicity in the absence of a systemic context is unknown. Here we demonstrate that while undifferentiated human ES cells do not possess an intrinsic functional clock, oscillatory expression of known core clock genes emerges spontaneously during directed cardiac differentiation. We identify a set of clock-controlled output genes that contain an oscillatory network of stress-related transcripts. Furthermore, we demonstrate that this network results in a time-dependent functional response to doxorubicin, a frequently used anti-cancer drug with known cardiotoxic side effects. Taken together, our data provide a framework from which the effect of oscillatory gene expression on cardiomyocyte physiology can be modeled in vitro, and demonstrate the influence of a functional clock on experimental outcome.

Cdx and T Brachyury Co-activate Growth Signaling in the Embryonic Axial Progenitor Niche.

In vertebrate embryos, anterior tissues are generated early, followed by the other axial structures that emerge sequentially from a posterior growth zone. The genetic network driving posterior axial elongation in mice, and its disturbance in mutants with posterior truncation, is not yet fully understood. Here, we show that the combined expression of Cdx2 and T Brachyury is essential to establish the core signature of posterior axial progenitors. Cdx2 and T Brachyury are required for extension of a similar trunk portion of the axis. Simultaneous loss of function of these two genes disrupts axial elongation to a much greater extent than each single mutation alone. We identify and validate common targets for Cdx2 and T Brachyury in vivo, including Wnt and Fgf pathway components active in the axial progenitor niche. Our data demonstrate that integration of the Cdx/Hox and T Brachyury transcriptional networks controls differential axial growth during vertebrate trunk elongation.

Polarized regulatory landscape and Wnt responsiveness underlie Hox activation in embryos.

Sequential 3'-to-5' activation of the Hox gene clusters in early embryos is a most fascinating issue in developmental biology. Neither the trigger nor the regulatory elements involved in the transcriptional initiation of the 3'-most Hox genes have been unraveled in any organism. We demonstrate that a series of enhancers, some of which are Wnt-dependent, is located within a HoxA 3' subtopologically associated domain (subTAD). This subTAD forms the structural basis for multiple layers of 3'-polarized features, including DNA accessibility and enhancer activation. Deletion of the cassette of Wnt-dependent enhancers proves its crucial role in initial transcription of HoxA at the 3' side of the cluster.

Transcriptional Dynamics at Brain Enhancers: from Functional Specialization to Neurodegeneration.

Over the last decade, the noncoding part of the genome has been shown to harbour thousands of cis-regulatory elements, such as enhancers, that activate well-defined gene expression programs. Driven by the development of numerous techniques, many of these elements are now identified in multiple tissues and cell types, and their characteristics as well as importance in development and disease are becoming increasingly clear. Here, we provide an overview of the insights that were gained from the analysis of noncoding gene regulatory elements in the brain and describe their potential contribution to cell type specialization, brain function and neurodegenerative disease.

Genome-wide DNA methylation levels and altered cortisol stress reactivity following childhood trauma in humans.

DNA methylation likely plays a role in the regulation of human stress reactivity. Here we show that in a genome-wide analysis of blood DNA methylation in 85 healthy individuals, a locus in the Kit ligand gene (KITLG; cg27512205) showed the strongest association with cortisol stress reactivity (P=5.8 × 10(-6)). Replication was obtained in two independent samples using either blood (N=45, P=0.001) or buccal cells (N=255, P=0.004). KITLG methylation strongly mediates the relationship between childhood trauma and cortisol stress reactivity in the discovery sample (32% mediation). Its genomic location, a CpG island shore within an H3K27ac enhancer mark, and the correlation between methylation in the blood and prefrontal cortex provide further evidence that KITLG methylation is functionally relevant for the programming of stress reactivity in the human brain. Our results extend preclinical evidence for epigenetic regulation of stress reactivity to humans and provide leads to enhance our understanding of the neurobiological pathways underlying stress vulnerability.

Epigenomic annotation of gene regulatory alterations during evolution of the primate brain.

Although genome sequencing has identified numerous noncoding alterations between primate species, which of those are regulatory and potentially relevant to the evolution of the human brain is unclear. Here we annotated cis-regulatory elements (CREs) in the human, rhesus macaque and chimpanzee genomes using chromatin immunoprecipitation followed by sequencing (ChIP-seq) in different anatomical regions of the adult brain. We found high similarity in the genomic positioning of rhesus macaque and human CREs, suggesting that the majority of these elements were already present in a common ancestor 25 million years ago. Most of the observed regulatory changes between humans and rhesus macaques occurred before the ancestral separation of humans and chimpanzees, leaving a modest set of regulatory elements with predicted human specificity. Our data refine previous predictions and hypotheses on the consequences of genomic changes between primate species and allow the identification of regulatory alterations relevant to the evolution of the brain.

Local compartment changes and regulatory landscape alterations in histone H1-depleted cells.

BACKGROUND: Linker histone H1 is a core chromatin component that binds to nucleosome core particles and the linker DNA between nucleosomes. It has been implicated in chromatin compaction and gene regulation and is anticipated to play a role in higher-order genome structure. Here we have used a combination of genome-wide approaches including DNA methylation, histone modification and DNase I hypersensitivity profiling as well as Hi-C to investigate the impact of reduced cellular levels of histone H1 in embryonic stem cells on chromatin folding and function. RESULTS: We find that depletion of histone H1 changes the epigenetic signature of thousands of potential regulatory sites across the genome. Many of them show cooperative loss or gain of multiple chromatin marks. Epigenetic alterations cluster to gene-dense topologically associating domains (TADs) that already showed a high density of corresponding chromatin features. Genome organization at the three-dimensional level is largely intact, but we find changes in the structural segmentation of chromosomes specifically for the epigenetically most modified TADs. CONCLUSIONS: Our data show that cells require normal histone H1 levels to expose their proper regulatory landscape. Reducing the levels of histone H1 results in massive epigenetic changes and altered topological organization particularly at the most active chromosomal domains. Changes in TAD configuration coincide with epigenetic landscape changes but not with transcriptional output changes, supporting the emerging concept that transcriptional control and nuclear positioning of TADs are not causally related but independently controlled by the locally associated trans-acting factors.

Inhibition of Super-Enhancer Activity in Autoinflammatory Site-Derived T Cells Reduces Disease-Associated Gene Expression.

The underlying molecular mechanisms for many autoimmune diseases are poorly understood. Juvenile idiopathic arthritis (JIA) is an exceptionally well-suited model for studying autoimmune diseases due to its early onset and the possibility to analyze cells derived from the site of inflammation. Epigenetic profiling, utilizing primary JIA patient-derived cells, can contribute to the understanding of autoimmune diseases. With H3K27ac chromatin immunoprecipitation, we identified a disease-specific, inflammation-associated, typical enhancer and super-enhancer signature in JIA patient synovial-fluid-derived CD4(+) memory/effector T cells. RNA sequencing of autoinflammatory site-derived patient T cells revealed that BET inhibition, utilizing JQ1, inhibited immune-related super-enhancers and preferentially reduced disease-associated gene expression, including cytokine-related processes. Altogether, these results demonstrate the potential use of enhancer profiling to identify disease mediators and provide evidence for BET inhibition as a possible therapeutic approach for the treatment of autoimmune diseases.

Large-scale identification of coregulated enhancer networks in the adult human brain.

Understanding the complexity of the human brain and its functional diversity remain a major challenge. Distinct anatomical regions are involved in an array of processes, including organismal homeostasis, cognitive functions, and susceptibility to neurological pathologies, many of which define our species. Distal enhancers have emerged as key regulatory elements that acquire histone modifications in a cell- and species-specific manner, thus enforcing specific gene expression programs. Here, we survey the epigenomic landscape of promoters and cis-regulatory elements in 136 regions of the adult human brain. We identify a total of 83,553 promoter-distal H3K27ac-enriched regions showing global characteristics of brain enhancers. We use coregulation of enhancer elements across many distinct regions of the brain to uncover functionally distinct networks at high resolution and link these networks to specific neuroglial functions. Furthermore, we use these data to understand the relevance of noncoding genomic variations previously linked to Parkinson's disease incidence.

Histone acetylation in astrocytes suppresses GFAP and stimulates a reorganization of the intermediate filament network.

Glial fibrillary acidic protein (GFAP) is the main intermediate filament in astrocytes and is regulated by epigenetic mechanisms during development. We demonstrate that histone acetylation also controls GFAP expression in mature astrocytes. Inhibition of histone deacetylases (HDACs) with trichostatin A or sodium butyrate reduced GFAP expression in primary human astrocytes and astrocytoma cells. Because splicing occurs co-transcriptionally, we investigated whether histone acetylation changes the ratio between the canonical isoform GFAPα and the alternative GFAPδ splice variant. We observed that decreased transcription of GFAP enhanced alternative isoform expression, as HDAC inhibition increased the GFAPδ∶GFAPα ratio. Expression of GFAPδ was dependent on the presence and binding of splicing factors of the SR protein family. Inhibition of HDAC activity also resulted in aggregation of the GFAP network, reminiscent of our previous findings of a GFAPδ-induced network collapse. Taken together, our data demonstrate that HDAC inhibition results in changes in transcription, splicing and organization of GFAP. These data imply that a tight regulation of histone acetylation in astrocytes is essential, because dysregulation of gene expression causes the aggregation of GFAP, a hallmark of human diseases like Alexander's disease.